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Unraveling the Jahn- Teller effect 
in Mn doped GaN 

A. Stroppa and G. Kresse 
Faculty of Physics, University of Vienna, and Center for Computational Materials Science, 
Universitdt Wien, Sensengasse 8/12, A-1090 Wien, Austria 

We present an ab-initio study of the Mn substitution for Ga in GaN using the Heyd-Scuseria- 
Ernzerhof hybrid functional (HSE). Contrary to semi-local functionals, the majority Mn t2 manifold 
splits into an occupied doublet and an unoccupied singlet well above the Fermi-level resulting in 
an insulating groundstate, which is further stabilized by a sizeable Jahn- Teller distortion. The 
predictions are confirmed using GW calculations and are in agreement with experiment. A transition 
from a localized to a delocalized Mn hole state is predicted from GaN to GaAs. 



Semiconductor based spintronics aims to develop hy- 
brid devices that could perform all three operations, 
logic, communications and storage within the same ma- 
terials technology Dilute Magnetic Semiconductors 
(DMSs) represent the most promising materials, and un- 
doubtedly, transition metal doped III-V semiconductors 
are presently the workhorse for spintronics Q. 

Ab-initio simulations based on density functional the- 
ory have played an important role in investigating the 
physics of DMSs 0,0]. Nevertheless, the theoretical un- 
derstanding has been hindered by the well-known defi- 
ciencies of the spin-polarized local density approxima- 
tion (SLDA) and generalized gradient approximation 
(SGGA) to the exchange-correlation functional the 
non-locality of the screened exchange interaction is not 
taken into account and the electrostatic self-interaction 
is not entirely compensated. This lack of compensation 
causes fairly large errors for localized states, e.g. the Mn 
d states. It destabilizes the orbitals and decreases their 
binding energy, leading to an over-delocalization of the 
charge density Another closely related issue is that 
the Kohn-Sham gap is usually a factor 2-3 smaller than 
the fundamental gap of the solidd. Whenever the energy 
position of the defect level with respect to the Valence 
Band Maximum (VBM) is comparable with the Kohn- 
Sham gap, e.g. deep acceptor levels introduced by sub- 
stitutional Mn in GaN (MnGa)!3|j the calculation of the 
thermodynamic transition levels becomes difficult, since 
all predicted thermodynamic transition levels are strictly 
bound by the Kohn-Sham one electron gap. Although the 
underestimation of the one electron gap would even occur 
for the exact Kohn-Sham functional, discontinuities in 
the potential upon adding or removing electrons correct 
for this error [^@]. For approximate functionals, which 
lack any such discontinuities — this includes all available 
semi-local and hybrid functionals — agreement between 
the Kohn-Sham gap and experimental fundamental gap 
is a prerequisite for modelling thermodynamic transition 
levels and band structure related properties Q. 

Hybrid Hartree-Fock density functionals [§| overcome 
the two limitations discussed above to a large extend^]. 



brid functional[lO| to study the Mn impurity in a GaN 
semiconductor host. Extensive studies of the perfor- 
mance of the HSE functional in solid state systems can 
be found in Refs.[il[i3, Il3l . unequivocally showing that 
hybrid functionals outperform semi-local functionals for 
materials with band gaps. The HSE results are confirmed 
by GWq calculations, which are the benchmark method 
for the prediction of quasiparticle (QP) energies [l3|. We 
will show that the electronic properties are accurately de- 
scribed by both methods. The ^2 manifold is split into an 
occupied doublet and an unoccupied singlet giving rise to 
a symmetry-broken insulating ground state that naturally 
couples with the ionic lattice, distorting the otherwise 
ideal tetrahedral environment of the Mn ion (Jahn- Teller 
effect). The calculated thermodynamic transition level 
e(0/— ) is in good agreement with experiments. Remark- 
ably, most of these features are not captured by standard 
SLDA or SGGA, without introducing ad-hoc corrections, 
such as self-interaction corrections 15, [3| or LDA+U 
corrections [17|. Hybrid functionals have a single param- 
eter (non-local exchange) that is once and forever fixed 
to 1/4 on theoretical grounds [l8|. Furthermore, Stengel 
et al. have very recently highlighted some drawbacks of 
SIC-schemes, concluding that hybrid functionals repre- 
sent the most promising route to reduce self-interaction 
problems while preserving unitary invariance[l9j|. 

The calculations were performed using the projector 
augmented- wave (PAW) method |20i with the Perdew- 



Burke-Ernzerhof (PBE) GGA functional[2l|land Heyd- 
Scuseria-Ernzerhof (HSE) hybrid functional l3| recently 
implemented in the VASP code^2], following exactly the 
prescription given in Ref. (HSE06). The Ga 3d and 
Mn 3p electrons were considered as valence electrons. A 
soft nitrogen PAW potential was used and the energy cut- 
off was set to 280 eV. Supercells with 64 and 128 atoms 
were used with lattice constants fixed to the optimized 
HSE value for the bulk crystal. Brillouin-zone integration 
was carried out using 8x8x8 and 2x2x2 Monkhorst- 
Pack grids for bulk GaN and the supercells, respectively. 
We used the Van der Walle and Neugebauer approach [24| 
for calculating the transition level of Mnca using a 128 
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TABLE L Lattice constant a, bulk- modulus Bo, energy gap 
at r, L, X, dielectric constant Soo, valence band-width W, 
and the energy position of Ga d states determined using 
PBE, HSE and GWo- Spin-orbit splitting is not included 
(Ao ~ 0.02 eV). £00 was calculated including local field ef- 
fects and using the RPA (values in parenthesis). 
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were relaxed. Errors due to the electrostatic interactions 
were taken into account through the Madelung energy 
of point charges q in an effective medium with a static 
dielectric constant £00 — 5-1, ^Ei—q^a/2eooL, where L 
is the distance between the Mn and its periodic replica, 
and a is the Madelung constant [2^. Soo was evaluated 
according to Ref. [l^. For the GWq calculations, we de- 
termined the screening properties entering Wq using the 
random phase approximation (RPA) and PBE wavefunc- 
tions and eigenvalues, but iterated the eigenvalues in G 
until selfconsistency was reached. The initial wavefunc- 
tions and eigenvalues in G were determined using the 
HSE functional. For bulk materials, this procedure yields 
results that are within 5% of experiment and the self- 
consistent scGW procedure with vertex corrections that 
was used in Ref. [2^. The latter procedure is out of 
reach for the systems considered here, but we expect the 
compromise to be very accurate, since PBE yields good 
screening properties in the RPA, and HSE excellent ap- 
proximations for the true QP wavefunctions [26j| . 

In Tab. [J we report the relevant equilibrium bulk pa- 
rameters of GaN in the zinc-blende phase. The HSE 
lattice constant is very close to the experimental value 
at T=0° K. The Ef^^ band-gap is 3.06 eV, almost twice 
as large as the PBE one (1.58 eV), again in very good 
agreement with experiment. Furthermore all HSE one- 
electron gaps are very close to the GWq QP energies. A 
slight discrepancy arises for the HSE valence band-width, 
for which PBE gives the best estimate (but experimental 
errors might be large). The position of the Ga d levels 
are described well using HSE, only surpassed by GWq. 
Finally, the HSE dielectric constant agree well with ex- 
periment. Overall, the agreement between experiment, 
HSE and GWq is very good. 

Let us now consider the Mn substitution at the Ga host 
site. Fig. [Ushows the band structure for the 64-atom cell 
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FIG. 1: (Color online) Spin-polarized projected band struc- 
ture of the Mn impurity in a cubic 64-atom cell using PBE 
(left) and HSE (right). Majority (minority) bands are plot- 
ted in the upper (lower) part. The horizontal line corresponds 
to the Fermi level (energy zero is equivalent to valence band 
maximum) . The blue (dark gray) and orange (light gray) cir- 
cles indicate the strength of the t2 (i.e. d^y+d^z+dyz) and e 
(i.e. d^2_y2 -1-^3^2 _j.2) character, respectively. The VBM and 
CBM of the host crystal are shown by small empty circles. 
The majority doublet and singlet QP-shifts are shown in the 
right part by blue triangles. 



lines A and S. In the PBE band-structure, the major- 
ity t2 states form an essentially dispersionless impurity 
band pinned at the Fermi level. They are three-fold de- 
generate at F with the partial occupancy of each state 
equal to 2/3. The three bands have predominantly t2 
character at the zone center (^^0.39), but they also show 
some anti-bonding contributions from the N-p states at 
the nearest N neighbors. The non-bonding e states can 
be found at the valence band maximum and they are 
strongly localized in the Mn sphere, with a total e char- 
acter of ^0.55. Below the e states, the GaN valence band 
states are found. They are slightly hybridized with Mn 
states. For the minority component, the t2 and e states 
are shifted above the host CBM due to exchange split- 
ting. The e states form a flat band even more localized 
than their majority counterpart (e character ~0.80). The 
Fermi level is located in the gap in the minority compo- 
nent, while it cuts the ^2 bands in the majority states. 
Therefore, in the PBE description, Mn2;Gai_a;N is a half- 
metal for a Mn concentration of x ^3.1 %. 

The HSE results differ significantly from the PBE re- 
sults. First, the band gap opening of the host is recog- 
nized. However, most relevant is that the majority ^2 



pied doublet, which remains above the host VBM. The 
energy separation of the singlet and doublet is 1.46 eV at 
r, and the system is clearly insulating. The singlet state 
is strongly localized in the Mn sphere (t2~0.40), whereas 
the doublet, which is close in energy to the valence band, 
is less well localized, with a d character of '^0.20. The e 
states do not form a flat band anymore, since they are 
pushed down in energy hybridizing more effectively with 
the host valence band. In the minority component, the 
t2 band is not as flat as in the PBE case (see Fig. [T]) 
but hybridizes more strongly with the first host conduc- 
tion bands, especially at the zone boundaries, whereas 
the e states remain mostly unchanged, apart from a rigid 
upwards shift. 

Remarkably, we observe that the quasiparticle GWq 
band structure is essentially identical to the HSE one- 
electron energies. In the GWq calculations, the doublet 
and the singlet shift upwards away from the VBM by 
0.25 eV and 0.43 eV compared to the HSE case, but the 
energy separation remains almost unchanged. This sug- 
gests that HSE is a legitimate shortcut for sophisticated 
many electron calculations, an observation that already 
transpires from the very good HSE one-electron band- 
gaps. The GWq calculations also clearly confirm that 
the metallic state observed in DFT-PBE is an artifact of 
the involved approximations. In fact, a band gap can be 
predicted using DFT-PBE as well, if (and only if) the 
band gap is calculated as the energy difference between 
the ionization potential and electron affinity calculated 
by removing and adding one electron in one supercell (~ 
1 eV). However, if a single electron or hole were placed in 
a huge supercell with many Mn substitutional sites, DFT- 
PBE predicts the electron affinity and ionization poten- 
tial to be equal (metal), whereas the hybrid functional 
predicts the band gap even in the limit of low electron or 
hole concentrations. The latter result is correct, whereas 
the PBE result is not in agreement with experiment. 

Another visible consequence of the splitting of the t2 
states is the Jahn- Teller distortion around the defect. 
The Mn t2 states do not transform as an irreducible rep- 
resentation of the Td group, since they are not 3-fold 
degenerate. This effect is usually termed "spatial symme- 
try breaking" [s^ . The splitting is even obtained at the 
ideal ionic structure, and for the ideal structure, three 
equivalent electronic solutions with identical energies are 
found. Each of the electronic solutions corresponds to a 
different charge ordering, and once the lattice is allowed 
to relax the nuclear framework distorts accordingly, i.e. 
the Jahn- Teller effect is at play. 

Fig. [2] shows the spin-density around the Mn atom in 
the (111) plane in the ideal structure, for both PBE (a) 
and HSE (b) calculated for a 128-atom cell, sampling the 
Brillouin zone at the F point only. The PBE spin-density 
is clearly symmetric, while in the HSE case the symmetry 
is spontaneously broken with the C3 symmetry around 
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FIG. 2: (Color online) PBE (a) and HSE (b) spin-density 
contours between -0.03 and 0.02 /^s/A'' plotted on the (111) 
plane. Blue (dark gray) large spheres are Mn atoms; green 
(light gray) large spheres are Ga atoms; orange (light gray) 
small spheres are N atoms. Red (dark gray) and blue (light 
gray) regions correspond to positive and negative spin-density. 
Note the lack of the C3 symmetry in (b) . 



group associated with the Jahn- Teller distortion, by con- 
sidering the character tables of all the possible subgroups 
of Tc;[34| and taking into account that i) the t2 manifold 
is split into a singlet and a doublet, and ii) that the sub- 
group must not contain the C3 symmetry operation. The 
only subgroup compatible with i) and ii) is T)2d- By al- 
lowing the ions to relax without symmetry constraints, it 
is found that the tetrahedral environment around the Mn 
atom becomes distorted: two of the Mn-N bonds relax 
to 1.98 A, while the other two bonds relax to 1.97 A (the 
ideal GaN distance is 1.94 A). The relaxed structure is 
within numerical uncertainty indeed consistent with the 
D2d symmetry, in agreement with experiments 1351] . The 
energy gain due to relaxations is 184 meV/Mn. Semi- 
local functionals (PBE) are not able to capture the Jahn- 
Teller distortion unless one introduces ad-hoc symmetry 
breaking displacements, but even then the energy gain 
due to the symmetry breaking is certainly significantly 
underestimated [36] . 

As a final confirmation of the HSE results, we have cal- 
culated the e(0/— ) thermodynamic transition level for 
MuGa in GaN by adding an electron in the HSE cal- 
culations and relaxing the geometry. We obtained a 
value of 1.9 eV comparing notably better with experi- 
ment (1.8±0.2 eV [33]) than PBE calculations (1.6 eV). 
In order to test the predictive capability of the pro- 
cedure, we have repeated the calculations for Mnca 
in GaAs, for which Mn is experimentally found in a 
d^ configuration[36]. In this case, the splitting of the 
t2 states drops to ~ 0.2 eV, and no Jahn- Teller ef- 
fect is theoretically found using HSE or experimentally 
expected [36j. Furthermore, the HSE band structure for 
Mn in GaAs (Fig. [3|) shows that HSE predicts a metal- 
lic behaviour for the majority electrons. Hence the hole 
introduced by Mn shows a transition from a localized to 
a delocalized character from GaN to GaAs, in full accor- 
dance with all experimentally available data. 

A further confirmation of the accuracy of the present 
approach is that it is capable to predict the position of 
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FIG. 3: (Color online) HSE spin-polarized projected band 
structure of the Mn impurity in a cubic GaAs 64-atom cell. 
See caption Fig[l]for details. 



ing energy of 4.2 eV in photoemission experiments for 
GaAsi38j, in very good agreement with the HSE results, 
while in SLDA or SGGA, the spectral weight is shifted 
towards the Fermi level. 

In summary, we have shown that the HSE functional 
predicts the localization of one electron hole at the Mn 
impurity in GaN giving rise to a Jahn- Teller distortion 
and an insulating ground state via splitting of the t2 Mn 
manifold. The existence of a band gap was confirmed by 
GWq quasiparticle calculations and is in agreement with 
experiment. The formalism is able to predict the relevant 
electronic features without imposing ad-hoc corrections 
on the d states and includes ionic relaxations. Further- 
more, a change of the electronic hole state from a local- 
ized character in GaN to an itinerant band in GaAs is 
predicted. This observation agrees with the fact that the 
localization is strongly environment dependent [l^ . Pre- 
vious interpretations of optical experiments [soj based on 
standard DFT electronic structure calculations should be 
revised according to the picture emerging in this study. 
Likewise, magnetic interactions should be re-investigated 
using functionals that predict a band gap in the quasipar- 
ticle spectrum. Extension to IV-group DMS are currently 
in progress, (ioj 
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